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ABSTRACT 

We consider implications of our new model of quasar lifetimes and light curves for the quasar luminosity 
function (LF) at different frequencies and redshifts. In our picture, quasars evolve rapidly and the lifetime 
depends on both their instantaneous and peak luminosities. The bright end of the LF traces the peak intrinsic 
quasar activity, but the faint end consists of quasars which are either undergoing exponential growth to much 
larger masses and luminosities, or are in sub-Eddington quiescent states going into or coming out of a period 
of peak activity. The "break" in the observed LF corresponds directly to the maximum in the intrinsic dis- 
tribution of peak luminosities, which falls off at both brighter and fainter luminosities. We study this model 
using simulations of galaxy mergers which successfully reproduce a wide range of observed quasar phenom- 
ena, including the observed column density distribution. By combining quasar lifetimes and the distribution 
of maximum quasar luminosities determined from the observed hard X-ray LF with the corresponding lumi- 
nosity and host-system dependent column densities, we produce the expected soft X-ray and B-band LFs. Our 
predictions agree exceptionally well with the observed LFs at all observed luminosities, over the redshift range 
considered (z < 1), without invoking any ad hoc assumptions about an obscured population of sources. Our 
results also suggest that observed correlations in hard X-ray samples between the obscured fraction of quasars 
and luminosity can be explained in the context of our model by the expulsion of surrounding gas due to heating 
from accretion feedback energy as a quasar nears its peak luminosity and final black hole mass. 
Subject headings: quasars: general — galaxies: nuclei — galaxies: active — galaxies: evolution — cosmology: 
theory 



1. INTRODUCTION 

The nature and evolution of the luminosity function (LF) 
of quasars at different redshifts and frequ encies has been 
studie d for more than thirty years (e.g., IScJunidj [19681; 
Schmidt & Gre enll983tlBovle et alJ2000UMivaii et all2000t 
Ueda et al. 20031 and references therein), but its relation- 
ship to the intrinsic properties of individual quasars is not 
well-understood. S pectral synthesi s modeling of the X- 
ray b ackground (e.g., iComastri et aljfl995t iGilli et aDll999l 
1200 11) as well as observed differences be tween hard X-ray 
and soft X-ray or optical quasar LFs (e .g., [Boyle et7i1fl99a 
lLa Franca et all2002UUeda et alJ20 03) imply (and require) a 
large population of optically obscured quasars. Ho wever, uni- 
fied m odels of active galactic nuclei (AGN) (e.g., Antonucci 
119931) which invoke geometric forms of obscuration as the 
dominant source of absorption cannot predict the distribu- 
tion of column densities or differences between LFs, but 
rather depend on these observations to determine the mod- 
eled form of obscuration. Even when calibrated by observed 
ratios of obscured to unobscured AGN, such models cannot 
account for measured quasar lifetimes or the selection-effect 
dependent differences in observed LFs at different frequen- 
cies and redshifts. Furthermore, a growing body of obser- 
vations imply isotropic or evolution-dependent obscuration 
which cannot be explained by these simple models alone 
(e.g.. lBorosodll992t iKuraszkiewicz e7al]l2000t lTrarJ | 2003l: 
iPase et alJ 12004 iBarger et alJ 120051 lAlexander et all l2005t 
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iStevens et all2005l) . 

Previous efforts to interpret the quasar LF have re- 
lied on restrictive assumptions about lifetimes and light 
curves of quasars, supposing, for example, that quasars 
either have universal lifetimes or that they evolve expo- 
nentia lly with time. Semi-analytic al modeling of t he LF 
(e.g., Kauffmann & Haehneltl 12000b IVolonteri et ail 12003b 
Wvithe & Loeb 2003) has neglected the obscured quasar pop- 
ulation and generally focused on reproducing the observed 
optical or soft X-ray LF, which not only has a different shape 
but also under-predicts the total quasar population by an order 
of magnitude at most redshifts and luminosities. 

Recently, we have begun to explore the impact of black 
hole growt h on galaxy fo rmation, using simulations of galaxy 
mergers ( Springel et al. 2005b). Our models reproduce the 
observed correlation between black hole mass and galaxy 
velocity dispersion (the Mbh - & relation) (|DiMatteoetalJ 
2005),' and link the quasar p hase of galaxies ( Hopkins et alJ 
2005a b) to galaxy evolution (Springel et al. 2005a). Further- 
more, the simulations predict qualitatively different quasar 
light curves than hav e been adopted in previous work 
(Hopkins et al. 2005a b). In our picture, the peak, exponen- 
tial black hole growth is determined by the gas supply over 
timescales ~ 10 8 yr, during which the gas inflows powering 
accretion generate large obscuring column densities. The 
growth shuts down when significant gas is expelled as it is 
heated by feedback from black hole accretion, creating a 
window during which the AGN is observable as an optical 
quasar for a lifetime ~ 10 7 yr, in good agreement with ob- 
servations, and yielding a significant obscured quasar popula- 
tion {Hopkins et al. 2005a). Hopkins et al. (2005b) analyzed 
simulations over a range of galaxy masses and found that the 
quasar light curves and lifetimes are all qualitatively similar, 
with both the intrinsic and observed quasar lifetimes being 
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strongly decreasing functions of luminosity, with longer life- 
times at all luminosities for higher-mass (higher peak lumi- 
nosity) systems. Moreover, they found that the resulting dis- 
tribution of column densities Nh depends significantly on the 
observed luminosity threshold, and agrees remarkably well 
with observed Nh distributions of both optical and X-ray sam- 
ples once the appropriate se lection effects are applied. 

In lHopkins et al we discuss the intrinsic distribu- 

tion of source properties obtained by applying our model to 
the quasar LF, recognizing the essential and realistic property 
that the time spent at a given luminosity depends on both that 
luminosity and the peak luminosity of the quasar (or, equiv- 
alently, the final black hole mass or host system properties). 
This results in a qualitatively different distribution of source 
properties than that implied by the idealized light curves that 
have been used earlier. 

Here, we combine our model of quasar lifetimes and the 
resulting distribution of intrinsic source properties with the 
luminosity and host system-dependent Nh distributions de- 
scribed above. With these self-consistent results derived from 
hydrodynamical simulations, we find that the typical column 
density distribution is a strong function of the instantaneous 
luminosity of a quasar, and fit it to simple analytical functions. 
Using the observed hard X-ray quasar LF to recover the dis- 
tribution of intrinsic source properties, we then combine our 
models of quasar lifetimes and the corresponding observed 
column density distributions to reproduce the expected LF at 
other frequencies given some absolute magnitude/luminosity 
limit. We find that our predictions for the optical B-band and 
soft X-ray LFs agree well with observations in both bands. 
Thus, our model of quasar evolution, without any assump- 
tions, naturally reproduces differences in the hard X-ray, B- 
band, and soft X-ray LFs over a range of redshifts. 

2. THE SIMULATIONS 

The simulatio ns presented her e are the series described in 
detailed in Hop kins et al J J2005bl) . performed with GADGET- 
2 (Springel 200 5), a new version of the parallel TreeSPH 
code GADGET (Springel. Yoshida. & White feOOll) based on 
a fully conservative formulation ( Springel & Hernauist 2002) 
of smoothed particle hydrodynamics (SPH), which is re- 
quired for energy and entropy to be simultaneously conserved 
when smoothing lengths evolve adaptively (see e.g., Hern- 
quist 1993, O'Shea et al. 2005). Our simulations account for 
radiative cooling, heating by a UV background (as in Katz et 
al. 1996b, Dave et al. 1999), and incorporate a sub-resolution 
model of a multiphase interstellar medium (ISM) to describe 
star fo rmation and supernova feedback ( Springel & Hernauist 
2003). Feedback from supernovae is captured in this sub- 
resolution model through an effective equation of state for 
star-forming gas, enabling us to stably evolve disks with arbi- 
trary gas fractions (see, e.g. Springel et al. 2005b; Robertson 
et al. 2004). 

Supermassive black holes (BHs) are represented compu- 
tationally by "sink" particles that accrete gas at a rate M 
estimated from the local gas density and sound speed us- 
ing an Eddington-limited prescription based on Bondi-Hoyle- 
Lyttleton theory. The bolometric luminosity of the black hole 
is Lboi = e r Mc 2 , where e r = 0.1 is the radiative efficiency. We 
assume that a small fraction (typical ss 5%) of Lboi couples dy- 
namically to the surrounding gas, and that this feedback is in- 
jected into the gas as thermal energ y. This fraction is a free pa- 
rameter, which we determine as in Di Matteo et al. (2005) by 
matching the observed Mbh - cr relation. For now, we do not 



resolve the small-scale dynamics of the gas in the immediate 
vicinity on the black hole, but assume that the time-averaged 
accretion rate can be estimated from the gas properties on the 
scale of our spatial resolution (< 30 pc). 

The progenitor galaxies in our merger simula- 
tions form a family with virial velocities V,,,> = 
80, 113, 160, 226, and 320kms _1 . The gas equation of 
st ate follows the multi-phase, star-forming structure derived 
in Springel & Hernauist (2003), resulting in a mass-weighted 
temperature of star forming gas ~ 10 5 K. For each simulation, 
we generate two stable, isolated d isk galaxie s , each with an 
extended dark matter halo with a iHernquisn dl99 0) profile, 
motivated by cosmological simulations (e.g. Navarro et al. 
1996; Busha et al. 2004) and observations of halo properties 
(e.g. Rines et al. 2000, 2002, 2003, 2004), an exponential 
disk of gas and stars, and a bulge. The self-similarity of 
any subset of these models is broken by the scale-dependent 
physics of cooling, star formation, and black hole accretion. 
The galaxies have masses M v ; r = V v 3 ir /(10GZ/o) with the 
baryonic disk having a mass fraction = 0.041, the bulge 
Wb = 0.0136, and the rest of the mass in dark matter with a 
concentration parameter 9.0. In Hopk ins et alJ l l2005al) . we 
describe our analysis of simulation A3, one of our set with 
V v tr = 160kms _1 , a fiducial choice with a rotation curve and 
mass similar to the Milky Way. We begin our simulation 
with pure gas disks, which may better correspond to the 
high-redshift galaxies in which most quasars are observed. 

Each galaxy is initially composed of 168000 dark matter 
halo particles, 8000 bulge particles, 24000 gas and 24000 stel- 
lar disk particles, and one BH particle. We vary the initial 
seed mass of the black hole to identify any systematic depen- 
dence of our results on this choice. In the cases considered, 
we choose the seed mass to be sufficiently small that its pres- 
ence will not have an immediate effect. Given the particle 
numbers employed, the dark matter, gas, and star particles are 
all of roughly equal mass, and central cusps in the dark matter 
and bulge profiles are reasonably well resolved (see Fig 2. in 
Springel et al. 2005b). The galaxies are then set to collide 
from a zero energy prograde orbit. 

3. COLUMN DENSITIES IN THE SIMULATIONS 
3.1. Determining Column Densities to the Quasar 

We determine the column de nsity between a black ho le and 
a distant observer as follows (Hopkins et al. 2005a b). We 
calculate the column density between a black hole and a hy- 
pothetical observer from simulation outputs spaced every 10 
Myr before and after the merger and every 5 Myr during the 
merger of each galaxy pair. We generate ~ 1000 radial lines- 
of-sight (rays), each with its origin at the black hole particle 
location and with directions uniformly spaced in solid angle 
d cos 9 d(j). For each ray, we then begin at the origin, calcu- 
late and record the local gas properties using GADGET, and 
then move a distance along the ray Ar = rjh sm i, where n < 1 
and /i sm i is the local SPH smoothing length. The process is 
repeated until a ray is sufficiently far from its origin (> 100 
kpc). The gas properties along a given ray can then be in- 
tegrated to give the line-of-sight column density and mean 
metallicity. We test different values of 77 and find that gas 
properties along a ray converge rapidly and change smoothly 
for 77 = 0.5 and smaller. We similarly test different numbers 
of rays and find that the distribution of line-of-sight properties 
converges for > 100 rays. 

Given the local gas properties, we us e the GADGET m ulti- 
phase model of the ISM described in Springel & Hernauist! 
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(2003) to calculate the local mass fraction in "hot" (dif- 
fuse) and "cold" (molecular and HI cloud core) phases of 
dense gas and, assuming pressure equilibrium between the 
two phases, we obtain the local density of the hot and cold 
phase gas and the corresponding volume filling factors. The 
values obtained co rrespond roughly to the fiducial values of 
iMcKee & Ostrikerl(ll977l) . A detailed description of the prop- 
erties of both cold and hot p hases of the resulting multi - 
phase structure can be found in Springel & Herna uistl ( 120031) . 
Briefly, cold phase (HI or molecular) clouds, with a tem- 
perature < 1000 K, contain a fraction > 0.9 of the total gas 
mass above a critical star-forming density threshold deter- 
mined by the equilibrium solutions to energy balance equa- 
tions for injection by supernova feedback and radiative cool- 
ing. Given a temperature for the warm, partially ionized com- 
ponent ~ 8000K (with densities n H ± ot ~ 0.1cm" 3 , or more 
accurately ~ 0.01 — 0.1n# total)* determined by pressure equi- 
librium, we further calculate the neutral fraction of this gas, 
typically ~ 0.3-0.5. We denote the neutral and total column 
densities as Nhi and Nh, respectively. Using only the hot- 
phase density allows us to place an effective lower limit on the 
column density along a particular line of sight, as it assumes 
a ray passes only through the diffuse ISM, with > 90% of the 
mass of the dense ISM concentrated in cold-phase "clumps." 
Given the small volume filling factor (< 0.01) and cross sec- 
tion of such clouds, we expect that the majority of sightlines 
will pass only through the "hot-phase" component. 

We assume t he in trinsic quasar continuum SED follows 
iMarconi et al.l (12004ft . based on optical through hard X- 
ray observations ("e.g.. lElvis et aljfl994t iG eorge et al1|1998[ 
I Vanden Berk et all2001tlPerola et al.|2002UTelfer et alJ2002t 
llJeda et alJl2003l IVignali et alJl2003l) . For the extinction at 
different frequencies, we consider a gas-to-dust ratio equal 
to that of the Milky Way, (A b /N H i)mw = 8.47 x 10 _22 cm 2 , 
but scaled by metallicity, A B /N H i = (Z/0.02)(A«/Nhi) -mw, 
as suggest ed by observations (e.g., iBou chet et al.l l 19851) . al- 
though Hopkins et al. (2005a) note that the resulting differ- 
ence is small. We use the Small Magellanic Cloud (SMC)- 
like reddening curve o fpeil dT992h . again motivated by obser- 
vations (Hopki ns et al.ll2004l) . We calculate extinction in X- 
ray frequencies (0.03-1 keV) using the photoelectric a bsorp- 
tion cross sections of Morrison & McCammon (1983) and 
non-relativistic Compton scattering cross sections, similarly 
scaled by metallicity. In estimating the column density for 
photoelectric X-ray absorption, we ignore the calculated ion- 
ized fraction of the gas, as it is expected that the inner-shell 
electrons which dominate the photoelectric absorption edges 
will be unaffected in the temperature ranges of interest. We 
do not perform a full radiative transfer calculation, and there- 
fore do not model scattering or re-processing of radiation by 
dust in the infrared. 

3.2. Column Density Evolution During Black Hole Growth 

For each simulation, we consider Nh values at all times 
with bolometric luminosity L = e r Mc 2 in some logarithmic 
interval, weighted by the total time along all sightlines a 
given Nh is observed. This then gives us a binned distribu- 
tion P(N H \L)dlog(N H ) cx tQogtNg), t(log(N H ) + dlog(N H ))). 
This distribution as a function of (observed) luminosity is 
shown in detail for our fiducial Milky Way-like simu lation 
(V vir = 160kms _I ) in Figure 3 of lHopkins et alJ ll2005bl> . in 
which we demonstrate that the resulting column density dis- 
tribution reproduces both the typical column density distribu- 
tion of hard X-ray selected quasar samples (e.g. JUeda etal.1 
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FIG. 1. — Lognormal median column density Nh and dis- 
persion ctn h as a function of luminosity and the ratio of in- 
stantane ous to peak luminosity, for the simulations described in 
Hopkins et al. 1 2005b). Points correspond to quasars with final black 
hole masses of 7 X 10 6 (black crosses), 3 X 10 7 (blue circles), 1 X 
10 8 (red squares), 3 X 10 8 (cyan stars), 7 X 10 s (green diamonds), and 2 X 
10'MQ(yellow triangles). The black line is the best-fit power law to Nh(L). 



120031) . and that of optically selected samples. In particular, 
the agreement with the column density distribution from the 
SDSS as determined in Hopkins et al. (2004) is very good, 
once similar selection effects are applied. 

At each intrinsic (un-attenuated) bolometric luminosity L 
we make a simple approximation to the observed distribution 
and fit it to a lognormal form, 

P(N H ) = 1 -= exp[-(log(Afe /Nh) 2 ) /(2<j 2 Nh )]. (1) 

CT« H v27r 

We show the resulting total (neutral plus ionized) Nh and 
<7jv w for all simulations as a function of L in Figure [2 Al- 
though we have shown this for the similar set of simulations 
described above, we have tested the column density distribu- 
tion as a function of luminosity across a wide range of simula- 
tions, varying the gas fractions, orbital parameters, gas equa- 
tions of state, concentrations, presence or absence of bulges, 
and seed black hole masses, and find that, although the fi- 
nal black hole mass (and corresponding peak luminosity) in 
any of these cases can be dramatically changed, the column 
density distribution as a function of instantaneous and peak 
luminosity shows no systematic dependence on any of these 
host properties (Hopkins et al., in preparation). Therefore we 
can adopt these fits with reasonable confidence across a wide 
range of redshifts and luminosities. Indeed, just the simulated 
luminosities above range from 10 7 - 10 14 Lq, covering the en- 
tire range of actual observed quasar luminosities at almost all 
redshifts. 

We find that the dependence of a^ H on L is weak, and 
we consider both constant o~ n h = 0.4 and a linear best-fit 
ctjVh = 0.7 + 0.1 \og{L/L pea k). It is important to note, however, 
that this gives only the dispersion for an individual simula- 
tion; in calculating statistical quantities such as the luminosity 
function, the dispersion across a population is needed. This is 
easy to determine based on the dispersion in Nh across simu- 
lations. Since we have assumed the distributions are individu- 
ally lognormal, the dispersion of the quasar population is sim- 
ply broadened to <r Ntl w 0.8 (or ct Nh = 1.0 + 0.1 log (L/L pefl ^)). 
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There is a clear trend of increasing Nh with L, which we fit to 
a power-law, giving 
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iV fl Ri4.2x 10 21 cirT 2 (L/10 u L Q ) 
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We note that the median neutral column density, Nh i, follows 
a similar relation, with a typical ionized fraction 0.3 — 0.5 
(mean 0.35). This form can be understood roughly as fol- 
lows, in the context of buried quasar growth during times 
when the black hole mass is growing to its final mass, be- 
fore peak quasar stages. Consider the time-dependent mass 
M c within the merging core of radius R c (~ lOOpc), and 
assume that the black hole grow s such that Mbh ~ T)M C 
(rj ~ 0.005) (M agorrian et alll998l) . The total density is then 
p = M c R~ 3 f(r), where f(r) is a dimensionless profile of or- 
der unity, and the column density is Nhi = (f*M c ) / (^mnR 2 ), 
where /* is the product of the mean neutral, hot-phase frac- 
tion (< 0.01 in the most dense regions of the galaxy) and the 
integral of the density profile f(r) (~ 1), and /i is the mean 
molecular weight. This gives Nhi ~ 10 21 cm" 2 mg, where 
M c = m% x 10 8 Mq. The luminosity, from a Bondi accretion 
model, is L = e r Mc 2 = 4irae r (GMBH) 2 pc 2 / c 3 , or using the def- 
initions above, L = (4Trae r rj 2 fo)(c/c s ) 3 (G 2 M 3 )/(cR 3 ), where 
fo = f( r ) ~ 1 as appropriate for some "accretion radius" and c s 
is the (mass-weighted) sound speed (~ 30km/s in the approx- 
imately isothermal core). These values give L as lO n L0mg. 

Thus, we expect N H i ~ 10 21 cm" 2 (L/10 11 L o ) 1 ' /3 , with varia- 
tion of the above parameters as well as dynamical and feed- 
back effects generating the considerable scatter seen in this 
relationship. This relation is much shallower than the naively 
expected relation Nh oc L expected if Mbh is constant (L oc 
p oc Nh) or L oc Mbh always, and strongly contrasts with uni- 
fication models which predict static ob scuration or N h inde- 
pendent of L up to some threshold (e.g., Fabian 1999). 

This modeling naturally produces a population of heavily 
attenuated objects at large luminosities. Although our fitted 
hot-phase column densities do not reach extremely Compton- 
thick levels Nh > 10 25_26 cm~ 2 , the brightest objects shown 
in Figure ^ reach median column densities ~ 10 23 cm" 2 or 
larger. As the fitted distribution is a lognormal about this me- 
dian with a dispersion across the quasar population a^ H ~ 1, 
this implies a significant population of heavily attenuated 



(Nh ~ IQ Z - 10 24 cm" 2 ) objects at soft X-ray (0.5-2 keV) 
frequencies; with the number density of sources falling of 
with large r Nh- This agrees well w ith both direct observa- 
tions (Treister et al. 2004; Mainieri et al. 2005) as well as syn- 
thesis models of the X-ray backgrou nd (Madau et alJl!994t 
IComastri et al J 1 995t lOilli et all 19991 120011) . which require a 
population of such objects, suggesting that our model should 
in principle account for the X-ray background spectrum based 
on a given luminosity function, without adopting arbitrary dis- 
tributions of source obscuration or additiona l obscured popu - 
lations. Furthermore, as discussed in Hopkins et al. (2005b), 
extension of this distribution of Nh(L) to very bright quasars 
in unusually massive galaxies or quasars in higher-redshift, 
compact galaxies which we have not simulated may, during 
peak accretion periods, reach Compton-thick values (Nh > 
10 25 cm" 2 ) of the typical column density. More likely, as 
our model assumes ~ 90% of the mass of the densest gas 
is clumped into cold-phase molecular clouds, a small frac- 
tion of sightlines will pass through such clouds and encounter 
column densitie s similar to those s hown for the cold phase 
in Figure 2 of [Hopkins et al. (2005a) and Figures 4 and 5 of 
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FIG. 2. — The column density distribution at three observed hard X- 
ray (2-10 ke V) luminosities from th e Milky Way like (V wr = 160kms~') 
simulation of Hopkins et al. 1 2005a), with a peak hard X-ray luminosity 
Lx, peak = 5 X lO^ergs -1 . The luminosities are chosen to be well below 
the peak luminosity (upper left) and near peak, corresponding to the final 
stages of obscured black hole growth (upper right) and subsequent "blowout" 
phase (lower left) as feedback from accretion energy expels gas in the cen- 
tral regions of the galaxy. Lower right shows the fraction of sightlines with 
Nh > 10 22 cirT 2 from the six simulation outputs closest in time and lumi- 
nosity to Lx. peak (black diamonds), with horizontal errors the spacing in 
luminosity and vertical errors corresponding to th e range across different 
possible column density prescriptions compared in Hopkins et al. 1 2005b). 
Solid line is the best-fit obscured fraction as a function of luminosity from 
ISimpsori<2005t) . Circles show the observations from Ueda et al. 1 20Q3), with 
the dashed line the expectation from their fitted Nh function. 



allows a large concentration of mass in sub-resolution obscur- 
ing structures, such as an obscuring toroid on scales < 100 pc, 
although many of the phenomena such structures are invoked 
to explain can be accounted for through our model of time- 
dependent obscuration. Geometrical effects may therefore, 
however, become relevant for interpreting and explaining the 
Compton-thick source population. 

3.3. Column Densities in Final Growth Stages 

Although the physical motivation for such a dependence 
is intuitive, the trend of increasing column density with 
luminosity seems to run o pposite to that ob s erved in a 
number of X-ra y sam p les lj$ teffen et al.l 12003 1 lUeda et alJ 
1 20031 lHasingerl 12004 iGrimes. Rawlines. & Willottl 120041: 
Sazon ov & Revnivtsevl 120041: iBarger et alJ 120051 ISimpson 
I200T In these samples, it appears that the fraction of broad- 
line objects and the fraction below some moderate column 
density (N H = 10 22 
with our prediction, 
veals that our simulations can and do, in fact, reproduce this 
behavior. Although it is possible to extend our spectral mod- 
eling of the quasar and column density calculation to describe 
the complete stellar population of the galaxy, and so deter- 
mine more accurately when an X-ray selected quasar will 
show a typical broad-line spectrum, we defer this to a later 
paper as it requires calculating colors and attenuations of the 
stars at all times in our simulations. However, we can easily 
compare to the observations of e.g., Ueda et al. (2003), and 
consider the fraction of objects with < 10 22 cm" 2 (more ac- 
curately, the ratio of number with Nh < 10 22 cm" 2 to all with 
N H < 10 24 cm" 2 ) as a function of hard X-ray luminosity. 



m 2 ) increase with luminosity, at odds 
However, a more detailed inspection re- 



HopkinsetaD^005b), as large as Nh ~ 10 cm . This also We examine the Nh distribution of our fiducial, Milky Way 
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like (V v ir = 160 kms -1 ) simulation in detail in Figure[2] Over 
most of the simulation, we find the general trend shown in 
Figure [2 and discussed above. The upper left panel shows 
the binned column density distribution (arbitrary scale) for 
all times with an observed hard X-ray luminosity LJ bs (2- 
lOkeV) = 10 41 - 10 42 ergs -1 (we use these units for ease of 
comparison with the observations quoted above), well be- 
low the peak observed hard X-ray quasar luminosity w 5 x 
10 44 ergs _1 , at which point the evolution is following the nor- 
mal trend shown in Figure^ However, when the quasar nears 
its final, peak luminosity, there is a rapid "blowout" phase 
as thermal feedback from the growing accretion heats the 
surrounding gas, driving a strong wind and eventually cut- 
ting off the accretion process, leaving a remnant with a black 
hole satisfying the Mbh - a relation ( Di Matteo et a fl l2005l 
Hocki nTet al.1 1200 ^). We analyze th is "blowout" phase in 
some detail in Hopkins et al. (2005a), and find that it can 
be identified with the traditional bright optical quasar phase, 
as the final stage of black hole growth with a rapidly de- 
clining density (allowing the quasar to be observed in opti- 
cal samples), giving typical luminositi es, column densit ies, 
and lifetimes of optical quasars (Hopkins et al 112005a b) in 
good agreeme nt with observations (e.g.. lHopkins et al.ll2004l 
Martinil2004l) . 

If we consider these stages of quasar growth, then, near the 
peak luminosity, we find a very different trend. The upper 
right panel of Figure |2] shows the column density distribu- 
tion for times with L£ bs w 10 44 ergs _1 , primarily in the final 
stages of the rapid obscured black hole growth phase, just be- 
fore feedback from the accretion begins to expel the surround- 
ing gas. The typical column densities are large, > 10 22 cm~ 2 , 
as we expect from our modeling above, although not high 
enough to generate large extinction in hard X-rays. How- 
ever, just ~ 10 Myr later, the quasar has expelled a significant 
amount of gas and column densities fall rapidly. The lower 
left panel of the figure shows the column density distribution 
for times with a typical luminosity L^ bs « 10 ergs" 1 , essen- 
tially the very peak quasar luminosity, corresponding to the 
bright observable phase in which the quasar is driving a wind 
and expelling nearby gas. The typical column densities are 
lower by an order of magnitude, Nh ~ 1 2 1 - 1 22 cm" 2 . In the 
lower right panel of the figure, we plot the fraction of sight- 
lines with column densities above 10 22 cm" 2 (the "obscured 
fraction"), for the six simulation outputs closest (in both time 
and luminosity) to the peak luminosity of the quasar. The re- 
sults are shown (black diamonds) with the range in luminos- 
ity bet ween each simu l ation ou tput shown as horizontal error 
bars. In lHopkins et alJ (l2005bl) we discuss in detail the differ- 
ences in column densities that result from varying our method 
of calculation, and find that (after accounting for some hot 
phase-cold phase separation) adopting various extreme cases 
yields a factor ~ 2 difference in the calculated column densi- 
ties, so we plot the resulting differences in the obscured frac- 
tion from this un certainty as vertical error bars. The data from 
lUeda et ail fc003l) (circles) and their expectation from their 
fitted column density distribution (dashed lines) are shown for 
comparison, as is th e best-fit obscu red fraction as a function 
of luminosity from Simpson (2005) (solid line). 

The agreement between the observations and our result is 
encouraging, especially as our calculation considers only one 
simulation, and is not necessarily meant to reproduce the trend 
in observed quasar populations. However, we do find a sim- 
ilar trend in all of our simulations near the peak luminosity 



(see also Hop kins et al.ll2005albl) . The key point is that we 
find, near the peak luminosity of the quasar as feedback drives 
away gas and shuts down accretion processes, the typical col- 
umn densities fall rapidly with luminosity in a manner similar 
to that obs erved. In our model f or the luminosity function, 
proposed in Hopkins et al. (2005c), quasars below the "break" 
in the observed luminosity function are either growing effi- 
ciently in early stages of growth or in sub-Eddington phases 
coming into or our of their peak quasar activity. Around and 
above the break in the luminosity function, the observed lu- 
minosity function becomes dominated by sources at high Ed- 
dington ratio at or near their peak luminosities. Based on the 
above calculati on, as well as the description of the blowout 
phase from Hopkins et al. (2005a), we then immediately ex- 
pect what is observed, that in this range of luminosities, the 
fraction of objects observed with large column densities will 
rapidly decrease with luminosity as the observed sample is in- 
creasingly dominated by sources at their peak luminosities in 
this blowout phase. This also further emphasizes that the evo- 
lution of quasars dominates over static geometrical effects in 
determining the observed column density distribution at any 
given luminosity. 

4. QUASAR LIFETIMES & THE LUMINOSITY FUNCTION 

Hop kins et al.l (2005c) showed that a proper accounting of 
realistic quasar light curves results in luminosity-dependent 
quasar lifetimes. In this picture, quasar lifetimes are functions 
of both the instantaneous luminosity and the peak luminosity 
(i.e. final black hole mass or host galaxy properties) of the 
system. Given a quasar lifetime above some luminosity as a 
function of the peak luminosity of the quasar, t(L' > L^L pea ^), 
the quasar LF (in the absence of selection effects) is given by 

I- dt(L Lpeak) d lp 

J d\og{L) 

where n(L pea k) is the rate at which of sources in a given log- 
arithmic interval in L pea k are "born" (created or activated) 
per unit volume and $(L) is the number density of sources 
per logarithmic interval in L. This formulation implicitly ac- 
counts for the "duty cycle" (the fraction of active quasars at 
a given time), which is proportional to the lifetime at a given 
luminosity. At any redshift, h(L pea i) will be, in general, a 
complicated function of the distribution of galaxy properties, 
including merger rates, masses, and gas fractions. However, 
having determined the quasar lifetime from our simulations, 
we can use an observed luminosity function to de-convolve 
ri(L pea k)- Since we are only interested here in demonstrating 
that our modeling self-consistently reproduces the observed 
differences in luminosity functions, in conce rt with the in- 
terpreta tion of the luminosity function from Hopki nTet alJ 
(2005c), we adopt this semi-empirical approach in what fol- 
lows, and thus all quantities in the equations determining the 
intrinsic and observed luminosity functions are completely 
determined. 

Given a distribution of values and some minimum ob- 
served luminosity L'"'", the fraction f„i, s of quasars with a 
peak luminosity L pea k and instantaneous bolometric lumi- 
nosity L which lie above the luminosity threshold is given 
by the fraction of Nh values below a critical N'^ ax , where 
Z™» = f v L exp(-<7„N% ax ). Here, fJL) = L v /L is a bolometric 
correction and <x v is the cross-section at frequency v. Thus, 

^>,i,C' ! )=^ln(^^), (4) 
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and for the lognormal distribution above, 



fobs(i / ,L,Lp ea i c ,L L , ) — 2 



1 + erf I 



■\og{N^T/N H ) 



V2 



(5) 



This results in a LF (in terms of the bolometric luminosity) 



$(^L,Lr)=i f f ohs (v,L,L peak ,L™") 
** J 

x dt ^' n(L peak ) d \og(L peak ). 
dlog(L) 



(6) 



The important point to recognize in this equation is that / &, 
has a complicated dependence on both instantaneous and peak 
luminosity (which does not simply factor out of the above in- 
tegrals), and thus the equations above will have a non-trivial 
dependence on the distribution of quasar peak luminosities 
n{L pea k) and the quasar lifetime, which are very different in 
our analysis from what has generally been used in previous 
modeling attempts. 

We consider qu asar lifetimes determine d from the simula- 
tions described in Hop kins et all fc005albl) . The light curves 
in the mergers are complicated, generally having a period 
of early rapid accretion after "first passage" of the galax- 
ies, followed by an extended quiescent period, then a tran- 
sition to a peak, highly luminous quasar phase, and then a 
dimming as self-regulated mechanisms expel gas from the 
galaxy center after the black hole reaches a critical mass 
and shut down accretion (Di Matteo et al. 120051) . While com- 
plex, [Hx>j>kin^uin (2005b) find that the total quasar lifetime 
?q(Z/ > L) above a given luminosity L is well-approximated 
by a truncated power law for every simulation studied, with 
t Q {U > L) = t 9 (L/\0 9 L Q ) a , where t 9 = t Q {U > 1O 9 L ) ~ 
10 9 yr, over the rang e 10 9 L^ < L < L nm {, for a given quasar. 
Given then that Hopkins et alJ ( 1200 5b) find this normaliza- 
tion is approximately constant across simulations, the life- 
time in each simulation is then entirely determined by the 
power-law slope a. The value of a depends on the peak lu- 
minosity of the quasar (equivalently, the final black hole mass 
or host system properties), with more massive (higher peak 
luminosity) quasars yielding shallower power-law slopes as 
they spend more time at high luminosities and large Edding- 
ton ratios. Over a wid e range of L pe( & (from ~ 10 10 - 1O 14 L ), 
Hopki nset alJ ll2005bl) find a = ct(L pea k) is approximately lin- 
ear with \ogL pea k, a = cto + a' \ogL pea k with an upper limit 
a = -0.2. The time spent in any logarithmic luminosity inter- 
val in this range is then simply 

dt/d\og{L)=\a\t 9 (L/\(?L & T. (7) 

Hopki nset al.l (12005 cl) examine a series of possible restric- 
tions to this model, which change ot(L pea k) slightly but yield 
very similar behavior, and we have considered all the cases 
described therein and obtain identical results in every case us- 
ing luminosity-dependent quasar lifetimes of this form. We 
show results for the case in which a is determined from the 
entire duration of our simulations, without imposing any arbi- 
trary cutoffs, which gives ao = -0.95 and a' = 0.32 for L pea k 
in units of 10 12 L Q ,i.e. a = -O.95 + O.321og(L pOTt /lO 12 L ). 

Given this model of dt/d\og{L) as a function of L and 
L pea k, we can then fit to any $(L) to determin e n{L nP „k)- The 
resultin g distributions fi(L pea i), discussed in Hopki ns et al.l 
(2005c) are fundamentally different from the naive expecta- 
tion of previous analyses of the LF which relied on ideal- 
ized models of the quasar lifetime, either assuming quasars 



"turn on" at a fixed luminosity for some universal lifetime, 
{dt jd log(L) oc S(L-L<x)) or assuming a pure exponential light 
curve over some interval {dt / d\og(L) = constant). These pre- 
vious models for the quasar lifetime yield a direct relationship 
between an observed luminosity and peak luminosity (final 
black hole mass), and predict a distribution of peak luminosi- 
ties n{L pea k) with essentially identical shape to the observed 
LF. However, accounting for the luminosity dependence of 
quasar lifetimes based on our detailed modeling, we find that 
quasars spend significantly more time at low luminosities than 
near their peak, resulting in a very different n{L pea k) distribu- 
tion. In this modeling, n(L pea k) has the same shape as the 
observed LF above the "break" in the luminosity function, 
and these quasars are accreting at high efficiency, near their 
peak luminosities. Below the break luminosity, however, the 
ri{L pea k) distribution turns over, and the faint end of the lumi- 
nosity function is dominated by sources near the break lumi- 
nosity (the peak of the n{L pea i) distribution) accreting either 
at high efficiency but in early stages of merger and growing to 
much larger luminosities, or in sub-Eddington phases going 
into or out of peak quasar activity (in the final stages of the 
host galaxy merger). The evolution of the LF with redshift, 
then, directly relates to evolution in n(L pea k), with the char- 
acteristic peak luminosity of quasars (final black hole mass 
being built) increasing with redshift as the break luminos- 
ity shifts to larger values. It is of particular interest to de- 
termine if this modeling of quasar evolution and the conse- 
quent novel interpretation of the luminosity function produce 
differences in the luminosity function in different wavebands 
(as well as column density distributions) consistent with what 
is observed. Furthermore, when comparing different mod- 
els of quasar lifetimes and the luminosity function, even for 
n(L pea k) distributions chosen such that two different quasar 
lifetime models will produce an identical luminosity function 
in a given waveband, any dependence of the column density 
distribution on peak luminosity (i.e. host system properties) 
or Eddington ratio will result in a different prediction for the 
luminosity function at all other frequencies. 

We consider the observed luminosity functions in the 
hard X-ray (<S> HX ; 2-10 keV), soft X-ray ($ sx ; 0.5-2 keV) 
and o p tical B-band (<j> B ; A B = 4400A), fro m lUeda et alJ 
(2003): lMivaii et al.l (2000): B ovle et alJ fcOOCl) . respectively. 
We rescale all luminosity functions to the same f2 m = 0.3, 
fl\ = 0.7, H() = 70km s" 1 Mpc" 1 cosmology. In order to 
make a direct comparison between luminosity functions, 
we further rescale all luminosity functions in terms of the 
bolometric luminosit y, using the bolometric corrections of 
iMarconi et alJ ( 120041) . It is important to note that using the 
con stant, lumin osity-independent bolometric corrections of 
e.g. Elvis et al. ( 1994) instead results in a significantly steeper 
cutoff in the luminosity function at high bolometric lumi- 
nosities, as the bolometric luminosity inferred for the bright- 
est observed X-r ay quasars i s almo st an order of magnitude 
smaller using the Elvis et al. ( 1994) corrections. However, it 
has been well-established that the ratio of bolometric lumi- 
nosity to hard or soft X-ray luminosity increases with increas- 
ing luminosity (e.g.. IWilk es et alJl l994t |Qreen et alJ Il995t 
IVignali et alJ2003tlStrateva et all2005l) . and furth er the sam- 
ple quasars of lElvis et al.Hl 9941) are X-ray bright ( Elvis et all 
2002). Accounting for the luminosity dependence of the UV 
to X-ray flux ratio, ceox, gives ri se to most of this di fference. 
We adopt the form for aox from Vignali et al] ( 120031) . but our 
results are relatively insensitive to the different values found 
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FIG 3. — Upper panel s: Hard X-ray (thick: lUeda et al.|2Q03|). soft X-ray 
(thin: IMivaii et alj |200(f), and B-band (dot-dash: lBovl eet alJ|200j3) LFs at 
z - 1.0, reseated using the bolometric corrections from Marco ni et all 120041) 
which we adopt throughout (left) and the luminosity-independent average 
corrections from Elvis et al. 1 1994) (right). Lower panels: The ratio of soft 
X-ray (solid) and B-band (dot-dash) LFs to the hard X-ray LF above. LFs 
are plotted ov er the corresponding range of the observed luminosities in 
the Ueda et al. 1 2003 ) hard X-ray LF at this redshift. Vertical lines bracket 
the minimum and maximum bolometric luminosities for which observational 
data exist (according to the appropriate bolometric correction) in the soft X- 
ray (solid line, only a minimum shown) and B-band (dot-dashed lines). Out- 
side of these ranges the fitted luminosity functions are extrapolated beyond 
the data. 



in the literature. These differences in the resulting bolometric 
luminosity functions are illustrated in Figure[3] The upper left 
shows the luminosity functions $//x (thick), $>sx (thin), and 
4>b (dashed) at z = 1 . 0, with bolometric luminosities from the 
Marc onTet al.l d2004) corrections and corresponding densities 
rescaled according to 

d$> _ d<£> d\ogL v 

— . (o) 

dlogL d\ogL y dlogL 

This can be compared to the same luminosity fun ctions con- 
verted using the constant bolometric corrections of Elvis et al. 
i 19941) (upper right). In both cases the qualitative differ- 
ences are similar, with the hard X-ray luminosity function 
significantly above the soft X-ray or optical LF at low lu- 
minosities just below the break, and the ratio of hard X- 
ray to soft X-ray or optical LF decreasing at and just above 
the break. This is seen in the lower panels of the figure, 
where we plot the corresponding ratios &sx / &hx (solid) and 
®b/&hx (dot-dash). However, it is apparent that the in- 
ferred number density from the X-ray LFs decreases much 
more slowly with luminosity using the Marconi et al. 
corrections which account for the luminosity dependence of 
aox, resulting in a larger gap at a given luminosity between 
the hard and soft X-ray LFs (and the optical as well). It is 
also immediately clear in this plot that the constant bolomet- 
ric corrections of Elvis et al. ( 1994) cannot apply uniformly 
to all luminosities and redshifts, as this actually predicts a 
larger number of optically selected bright quasars than soft 
or hard X-ray objects, which cannot be explained with any 
sort of reddening/obscuration model. This is not a fitting 
artifact, as a direct com parison of the data between, e.g., 
ICroom et ail <2004l) and Barker et all J2005I) in the optical 
and h ard X-ray, respectively, shows (using the Elv is et aD 
i 19941) corrections) the optical quasar LF to be a factor ~ 2— 3 



higher than the hard X-ray LF at several luminosities and red- 
shifts, and this probl em is only w orse when consider i ng the 
optical LFs of IBovle et all (120001) or IRichards et all (120051) 
which are steeper at low luminosity than that of Croo m et alJ 

(2004) . Furthermore, there is direct evidence of a signifi- 
cant nu mber of optically obscured quasars, eve n at high lumi- 
nosities ( Norman et al. |200llStern et alJl2002t lDawson et alJ 
120031 iTreister et alJl2004l) . many of which are un dergoing a 
buried quasar or pre-quasar growth phase (iPa ge et alJl2004t 
lAlexander et all2005tlStevens et al.l2005l) as predicted in our 
modeling. As discussed above, the existence of a signifi- 
cant population of such objects is also inferred from syn- 
thesis models of the X-ray background (|Madau et alJ 119941 
IComastrietaill 19951 iGilli et alJI 19991 1200 ll) . At theextreme 
low and high luminosities shown, the relative behavior of the 
plotted LFs becomes confusing, but this is because they are 
extrapolated well beyond the range of observed data. There- 
fore, we show the actual range of observed luminosities for 
each LF as vertical lines (of the same style as the correspond- 
ing LF). The difference in LFs across these ranges are much 
less dramatic, and this is what we are interested in making 
a detailed comparison with. It is clear, though, that it is im- 
portant to account for the luminosity dependence of quasar 
bolometric corrections, as it creates this significant difference 
in the high-luminosity end of the bolometric quasar luminos- 
ity function and implies that a non-negligible fraction of the 
brightest quasars are not seen in optical surveys, and further 
that the ratio of one luminosity function to another is not triv- 
ially related to the obscured fraction discussed in detail above. 
For fu rther comparison of t hese bolometric corrections we re- 
fer to [Marconi et al. (2004), and for further detailed compar- 
ison of luminosity functions using a very si milar procedure 
to prod uce the bolometric corrections see, e.g. . IRichards et ail 

(2005) . 

Figure |4] shows the resulting LFs in different bands at red- 
shifts z = 0.3,0.5,0.8, and 1.0. All quantities have been 
rescaled in terms of the bolometric luminosity as discussed 
above. We calculate the n(L pea k) distribut ion by fitting to th e 
observed hard X-ray (2-10 keV) LF of lUeda et all ( I2003I) . 
<3?hx at each redshift. Given this n{L pea k), we then use our 
model of quasar lifetimes and the Nh distributions determined 
in §2 to calculate the expected B-band (<5> B ; X B = 4400A) 
and soft X-ray (&sx\ 0.5-2 keV) LFs, given the appropriate 
redshift-dependent sample luminosity/magnitude limit L"™ 
(simply taken as the minimum observed luminosity of each 
sample at that redshift). We comp are these pre dicted LFs to 
the corresponding observed IBovle et alJ J2000) B-band and 
Mivai i et alJ ( 12000) soft X -ray LFs at these redshifts. The 
agreement between the observed and predicted LFs is excel- 
lent at all observed luminositie s, and is repro duced for all low 
redshifts modeled. In Hopkins et al J ( 1200 5b), we have consid- 
ered several approaches to calculate the column density dis- 
tribution, and find that, after accounting for the clumping of 
most mass in some hot phase-cold phase separation, consid- 
ering or ignoring metallicity and ionization results in factor 
~ 2-3 differences in the calculated column densities and cor- 
responding quasar lifetimes along a given line of sight. Thus, 
we consider this to be a rough parameterization of the ex- 
tremes of our modeling, and estimate the resulting range in 
the predicted luminosity functions as a consequence of these 
extremes in calculations of the column density. We show the 
typical (averaged over the plotted points) resulting range in 
dex for both the optical and soft X-ray LFs at z = 1 .0 in the up- 
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Fig. 4. — Hard X-ray (thick), soft X-ray (thin), and B-band (dot-dash) 
LFs determined from our model of quasar lifetimes and column densities, 
based on a distribution of intri nsic source properties fitted to the observed 
hard X-ray LF as in Hopkins et al. 1 2005c) and the limiting magnitudes of 
observed samples, at the different redshifts shown. All quantities are rescaled 
to bolo metric luminosities with the bolometric corrections of Mar coni et all 
120041) . Symbols sho w the observed LF s over the range where observations 
exist, for hard X-rays i Ueda et al. 2003, diamonds), soft X-rays i Mivaii et al. 
2000 , triangles), and B-band i Boyle et al. 2000 , crosses). In the z = 1 panel, 
the points at upper right show the mean systematic offset in the soft X-ray 
(open) and B-band (filled) luminosity functions which result from modi- 
fying our column density prescription to ignore ionization, metallicity, or 
vary the attenuation calculation following the range of possibilities studied in 
Hopki ns et al] <2005bl) . 



per right of the corresponding panel, and note that the range at 
the other plotted redshifts is similar (though slightly smaller). 
The plotted range is significant in considering the relative lu- 
minosity functions, but we plot them not as errors but as an 
upper limit to the systematic effects of various extreme as- 
sumptions such as ignoring ionization and metallicity in cal- 
culating the column density and quasar attenuation. Based on 
our analysis above showing that near very peak luminosities 
the quasar will drive a wind expelling nearby gas and rapidly 
reducing the column density, we also consider a column den- 
sity distribution which follows the above, fitted form for most 
of the quasar lifetime but cuts off (obscuration is neglected) 
when the quasar is within a fraction ~ 1/2 of its peak lumi- 
nosity. We find that this makes little difference to the pre- 
dicted luminosity functions, primarily because the time spent 
at these very near peak luminosities is relatively small, and 
further because this results in a somewhat compensating small 
shift in the fitted n(L pea k) distribution. We note again though 
that a more complete treatment of this effect requires more 
complete spectral and dust modeling of both the quasar and 
stellar populations of the host galaxy, which we defer to a 
later paper. At higher redshift z > 1, we recognize that the 
light curves and Nh distributions may evolve as a result of 
changing host galaxy properties, and we defer a modeling of 
the LFs at high redshifts to a future paper. 

We also consider the results obtained using our column den- 
sity distributions, but applying only the idealized, luminosity- 
independent models of the quasar lifetime described above. 
The difference in lifetim es and the r esul ting h{L pea k) are de- 
scribed in detail in Hopkins et al. (2005c), but essentially 
n(Lp ea k) oc $(L = L pea k) for these models. For the simplest 
fit to the Nh distributions in §2, N H = N H {L) oc L° 35 and 
<tn h = constant « 0.8, f b s is independent of L pea k and can be 



taken out of the integral, giving $(i/, L, L™") = f obs (i/,L) <1>(L), 
independent of the lifetime and n(L pea k) model. However, 
even for a weak dependence on L pea k, o^ H = a^ H (L,L pea k) = 
1.0 + 0.1 log (L/ Lpeak), we find that using these models of the 
quasar lifetime under-predict both $b and $sx by a factor of 
> 3 at low and high luminosities. This is because these mod- 
els do not account for the quasar spending most of its life 
at luminosities well below its peak and thus do not properly 
account for quasars with different L pea k (i.e. different host 
galaxy properties such as total mass or gas fraction) at a given 
observed luminosity. In any case, such a procedure is not 
self-consistent, as the data from which our NH{L : Lp ea k) and 
aN H (L,Lp ea ic) relations are fitted imply and produce our model 
of luminosity-dependent quasar lifetimes, with the vast ma- 
jority of each Nh{L) distribution corresponding to points on 
the lightcurve which do not exist in these idealized models. 

5. CONCLUSIONS 

Using our picture of merger-driven quasar activity with 
self-regulated black hole growth and feedback, we are able 
to simultaneously reproduce the observed hard X-ray, soft 
X-ray, and B-band luminosity functions (LFs) over a broad 
range of observed luminosities and redshifts with significantly 
greater accuracy than previous models and without invoking 
any assumptions beyond the basic input physics of our sim- 
ulations. Furthe rmore, our picture als o yields the observed 
M BH -a relation (Di Mat teo et alJ2005l). the bimodal distribu- 
tion of galaxy colors (Srjrinael et al. 2005a), observed quasar 
lifetimes (Hopkins et al. 2005a), t he Nh dis tribution of both 
optical and X-ray samples (Hopkins et al. 2005b), and the 
faint-end slope of the quasar LF and sup ermassive black hole 
mass distribution (Hopkins et al. 2005c). We note that we are 
not predicting the luminosity function in this work, but are 
demonstrating that the observed differences between lumi- 
nosity functions are essentially entirely accounted for in our 
modeling. In other words, given a luminosity function in any 
waveband, our modeling allows us to accurately determine the 
luminosity function in other wavebands, giving us robust pre- 
dictive power across different samples and demonstrating that 
these observed differences can be explained self-consistently 
through a model in which observed patterns of quasar obscu- 
ration are dominated by differences in different stages of the 
evolution of quasars, and not simply by viewing-angle effects. 

Our model predicts that self-regulating feedback processes 
in galaxy mergers reproduce the difference in the quasar LF at 
different frequencies naturally, as a consequence of the evolu- 
tion of gas flows fueling accretion from gravitational torques 
(e.g. Barnes & Hernquist 1991, 1996; Mihos & Hernquist 
1996) and accretion feedback. The population of obscured 
quasars is also a natural consequence of the model, not as an 
independent population but as a stage in the "standard" evo- 
lution of quasars over their lives, before feedback can clear 
sufficient material to render the quasar visible. Once feed- 
back begins to unbind gas in the central regions of the galaxy, 
the quasar becomes observable in optical wavebands, and col- 
umn densities rapidly decrease near the quasar peak luminos- 
ity. This suggests that our modeling of this blowout phase, 
coupled with the new interpretation of the luminosity func- 
tion resulting from our model of quasar lifetimes, can explain 
the observed trends in the fraction of obscured or broad-line 
quasars with luminosity. 

The close agreement between our predictions and the ob- 
served LFs is strong evidence in favor of our self-consistent 
model of quasar lifetimes and light curves. This model sug- 
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gests a new and qualitatively different interpretation of the 
quasar LF, which we propose in Honkins et al. (2005c). Our 
interpretation of the quasar LF and the intrinsic deconvolved 
distribution of peak quasar luminosities and host galaxy prop- 
erties has important implications for the evolution of quasar 
populations, the energetics of the cosmic X-ray, UV, and IR 
backgrounds, the role played by quasars in reionization, and 
the production of the present-day distribution of supermas- 
sive black holes. Future attempts to understand, model, or in- 
corporate the distribution of quasar properties should account 
for the difference between the observed LF and the intrinsic 



distribution of source properties as a result of luminosity de- 
pendent quasar lifetimes, and the simultaneous, luminosity- 
dependent effects of evolving Nh distributions. 
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